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Abstract 

An ensemble of multilevel atoms is a good candidate for a quantum information storage device. 
The information is encrypted in the collective ground state atomic coherence, which, in the absence 
of external excitation, is decoupled from the vacuum and therefore decoherence free. However, 
in the process of manipulation of atoms with light pulses (writing, reading), one inadvertently 
introduces a coupling to the environment, i.e. a source of decoherence. The dissipation process is 
often treated as an independent process for each atom in the ensemble, an approach which fails at 
large atomic optical depths where cooperative effects must be taken into account. In this paper, 
the cooperative behavior of spin decoherence and population transfer for a system of two, driven 
multilevel-atoms is studied. Not surprisingly, an enhancement in the decoherence rate is found, 
when the atoms are separated by a distance that is small compared to an optical wavelength; 
however, it is found that this rate increases even further for somewhat larger separations for atoms 
aligned along the direction of the driving field's propagation vector. A treatment of the cooperative 
modification of optical pumping rates and an effect of polarization swapping between atoms is also 
discussed, lending additional insight into the origin of the collective decay. 
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I. INTRODUCTION 



A system of atoms interacting with a common reservoir (the electromagnetic vacuum) 
is often treated using an assumption of independent dissipation rates for the atoms. This 
is a valid assumption for the case when the interatomic distances are large (compared to 
an optical wavelength); however, when the distances between atoms become smaller or 
comparable to an optical wavelength, the mode structure around one atom is changed due 
to the presence of other atoms located in its immediate vicinity, and the decay rates are 
modified. Equivalently, the radiation emitted by one atom can be scattered off the second 
atom, thus changing the radiative properties of the system. As a consequence, the radiative 
decay of the ensemble must be viewed as a cooperative effect. A quantitative analysis of 
cooperative effects was given by Dicke [1] for an ensemble of two- level systems confined to 
a spherical volume whose radius is much smaller than an optical wavelength. For such a 
system, the collective decay rate can be increased by a factor proportional to the number of 
atoms in the ensemble. Further investigations extend the treatment to arbitrary interatomic 
distances, although the calculations become more complex. Formalisms for the treatment of 
cooperative spontaneous emission and resonance fluorescence from a system of many atoms 
have been developed [2-6] , and have been applied to a system of two two-level atoms [7-9] . 

Cooperative decay in multilevel atomic systems has yet to be treated in detail. A mul- 
tilevel atom has properties not possessed by a two-level atom: for example, it can store 
information in superpositions of ground state magnetic sublevels. Such ensembles are ex- 
tensively discussed in the literature as convenient systems for information storage or large 
scale entanglement generation [10]. In particular, pencil-shaped media have been used for 
the generation of spin squeezed and Schrodinger cat states, in the context of continuous 
measurement of a scattered field [11]. In such schemes, the coupling to the vacuum has a 
two-fold function: on one hand it gives rise to the signal while, on the other hand, it leads to 
an irreversible leakage of information from the system to the environment. In treating the 
losses due to spontaneous emission, the above mentioned assumption of independent atoms 
is generally used, which is a sound assumption as long as the atomic density is low. How- 
ever, optimal results (e.g. strong entanglement) are found in the regime of resonant optical 
depths greater than unity, a regime in which pencil-shaped media of two-level atoms exhibit 
superradiant behavior [12]. Even if the dynamics of a collection of multilevel atoms might be 
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substantially different from that of the two-level ensemble, the validity of the independent 
spontaneous emission regime is questionable. 

We proceed in the present publication with an analysis of cooperative effects in a system 
of two, four-level atoms. This provides the starting point for an extension to many atom 
systems, while it also addresses the non-trivial question of the importance of cooperative 
decoherence in a simple quantum information system of two qubits. The calculations are per- 
formed foraJ=l/2^J' = l/2 transition irradiated with a monochromatic, off- resonant, 
<7 + polarized, classical laser field. The collective decoherence of an initial equal superpo- 
sition of ground sublevels (x polarized atomic state) is obtained for arbitrary interatomic 
separations and compared with that of atoms independently coupled to the reservoir. In 
addition, the transfer of population from a z polarized atomic state with both atoms in one 
of the ground sublevels to another z polarized state with both atoms in the other sublevel, 
is analyzed. A polarization swap effect is also discussed, where an x polarized atom induces 
x coherence in a neighboring atom [13]. 

The paper is organized as follows: in Sec. II the theoretical method used is described. In 
Sec. Ill analytical solutions for the case of close atoms are obtained. Numerical solutions 
for arbitrary separations are discussed and plotted in Sec. IV. In Sec. V the polarization 
swap effect for arbitrary separations is discussed, while Sec. VI contains some conclusions. 

II. THEORY 

As indicated in Fig. 1(a), the two atoms (having natural frequency u ) are located at the 
origin, Ri = 0, and at R2 = R21, respectively. The traveling wave driving field propagates 
in the positive z direction with wave vector k, frequency Q (detuned from ujq by A), and 
circular polarization a + . Denoting the 4 states of a single atom with | {) (ground, m| = — 1/2 
eigenvalue), \]) (ground, = 1/2 eigenvalue), \(3) (excited, mp = —1/2 eigenvalue) and 
I a) (excited, m a =1/2 eigenvalue), the classical field drives the ||) — > \a) transitions. The 
Hamiltonian for the system (composed of the two atoms indexed with i, where i — 1, 2) is a 
sum of the free Hamiltonian (H ), the classical field-atom interaction Hamiltonian (V c ) and 
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FIG. 1: (a) The classical field is propagating along the z axis, driving the two atoms, located at the origin 
Rl = and at R2 = R21, respectively, (b) Internal structure of a single atom. 



quantized bath-atom interaction Hamiltonian (V^) given by: 



#0 = E ^0 [\oi)i (a\ { + 

i=l,2 



+ E^fc a kA a k,A 
k,A 



(1) 



V c = E [|a>< (II, e M e-^ + 1 1). H e~ M e* m ] 

2=1, r 

K = E 



8=1,2 



E E E* <& I^Oli <Ae-* k ' R * +« A ) I/}, 01, a M e' 



_i=i,Ti'=/3,ak,A 

As usual, the field is quantized in a volume V and is described by the creation and annihi- 
lation operators a} kX and a^\, corresponding to modes with wave vector k and polarization 
A. The atomic dipole moment d (the same for both atoms) couples to both the classical 
and quantized fields. The cw driving field at the position of the atoms (R, for £ = 1,2) 
is expressed as E (R, t) —\ [E e* kL Ri e~ %nt + cc] e+ . The classical part of the interaction 
contains the Rabi frequency defined as: 

d + E 



X 



2h 



(2) 
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where d + is the matrix element d + = (a\ d-fT+ ||); the interaction with the quantum vacuum 
has an associated coupling strength 

proportional to the dipole matrix element in the direction of the unit polarization vector 
<1 =0'|d-e M \j>). 

The quantities that are relevant in what follows are the collective coherence operator 

At = E U>< <TU (4) 

i=l,2 

and the normalized population operator 

At = ^E IT>* <TU • (5) 

The derivation of the time evolution of the expectation values of these two operators is the 
goal of our calculations. They can expressed in terms of density matrix elements as 

(At ) = PTT;TT +PTUT +PTT;TT + AlT;U > ( 6 ) 
(At ) = 2 ( 2 ^TT;TT + PfUU + PTT;TT ) ■ 

We consider two problems: first, one with both atoms prepared initially in a superposition 
of ground states with maximum coherence ^ + |T)) (collective coherence equal to 1) 
and, second, one where the population is transferred from the state with both atoms in the 
||) state to the one with both atoms in the |f) state. 

Qualitative (and some quantitative) details of the calculations are described in the follow- 
ing. The subspace of interest in which the collective operators defined above act (henceforth 
named the ground subspace) is of dimension 4 and it is spanned by state vectors containing 
the ground substates of the two atoms ( |||),||T), IT I) and |TT))- A set of 16 density matrix 
equations completely describes the dynamics of this space. However, the ground subspace is 
coupled to the ground-excited subspace of dimension 8 (containing states with one excitation 
as for example |J, a)) through the classical field. This is, in its turn, coupled to the excited 
subspace of dimension 4 (containing states of two excitations like \aa)) which can decay 
back to the ground states. In a density operator approach, a total of 256 density matrix 
elements coupled to each other come into play, which makes the task at hand extremely 
complex. 
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Some simplifications are possible. First, terms occurring in the evolution of the ground 
state density matrix elements are separated into in-terms (due to spontaneous emission 
from upper states) and out-terms (driving terms due to the presence of the classical field), 
and these terms are treated separately. Second, an amplitude rather than a density matrix 
approach is sufficient to obtain expressions for the excited-ground and excited subspace 
density matrix elements which enter the equations. The procedure is described in Appendix 
B where it is applied to the derivation of the decoherence and population transfer rate for a 
single 4-level atom. The treatment is perturbative in the sense that excited state populations 
are assumed to be negligibly small, as in many treatments of optical pumping. 

The states coupled by the fields in this approximation are denoted by \fjfu), \^v') and 
\/j,u) where the convention used is that the prime indicates excited states (a or f3), while 
the unprimed symbols indicate ground states (j or |). By eliminating the intermediate 
states involving the radiation field (procedure outlined in Refs. [14, 15]), one can obtain 
the following coupled equations of motion for the state amplitudes containing one excitation 
(b^i and v^)> in an interaction picture: 

(R21) {m M , Try } { m » » «v} (7) 

-iAb^ +i X e ikL - ni b iv 5 m ^,, l/2 , 

bf,v> = -7V' ~ lJ2 G m„ -m^,;m v -«v (R12) {"V > "V } {"V ' m ^ } V" 

- iAb^ + ?xe iki ' R2 b^ <5 lA nv • 

In the equation above for b^ u , the first term on the right-hand side is the decay (at a 
rate 7 equal to half the excited state population decay rate) of the excited state amplitude 
of an atom independently coupled to the quantum vacuum. The second term contains 
a propagator G m ^- m ^ m ^- m „ (R21) which includes the effects of the radiation exchange 
between atoms: an atom makes a transition from state v' (quantum number m v i) to state v 
which is accompanied by a transition in the other atom (// to //). The real part of G gives a 
contribution to the decay rate that varies from 1, for maximum cooperation between atoms 
when their separation is much less than A, to 0, when no exchange of radiation between 
atoms is present (infinite separation). The imaginary part leads to a shift in energy (which 
adds to A in the equations above) and varies from (large separation) to infinity when the 
atoms are in the same location. If the minimum interatomic separation is small but finite, 
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the shift can always be kept small compared to the detuning and can be neglected. The 
explicit expressions for the propagators involved in this problem are given in Appendix A. 
The geometrical information on the radiation exchange is contained in the Clebsch-Gordan 
coefficients {m,m'} = (1/2, 1/2; m, m'\l, m — m'). The last two terms in the right-hand 
side arise from driving field induced transitions and from the off-resonant nature of the 
interaction. 

To calculate the out-terms, we note that, as a result of the driving field, ground state 
amplitudes {b^ v ) are coupled to excited state amplitudes via 

K v = ix [e^ k - Rl h av <L lAmM + e- ik ^ R2 b, a <L 1/2 , m „ ] . (8) 

As a result one finds 

Pav;mn = Knn + ^ ^mn = [^X e % L * ^au 5_l/2,m„ + e 2 £-l/2,m„ J Knn + 

(9) 

+ [i X e- lkL ' Rl b an 5_ 1/2>rnm + e - lkL - R2 b rna 5_ 1/2)mn ] . 

The system of equations (Eqs. (7)) is solved for the 8 amplitudes b^ v and b^ v > as functions 
of the 4 ground state amplitudes b^ u ; these expressions are replaced in the above equation 

. out 

and with the identification b^ v b* mn — > p^ mn , the rate equations for p^ v . mn are obtained in 
terms of the 16 density matrix elements p M ^ ;m „ (with /i, u,m,n —I, |). 

Next, repopulation from the upper states to the lower states is taken into account (in- 
terms) 

Pmnifiu = ^m/i Pm'n;fi'u + X] ^nu Pmn';^u' + (10) 

m' ,n' n' ,u' 

( 2 7) Yl G mi// - mv . imml - mrn (R 2 i) {m m ,m m > } {m u ,m v > } p m 'n^v' + 

m' ,v' 

(R21) {m M , Try } {m n , m n >} 

(11) 

m' ,u' 

where 

r ab 6 ' = ( 2 7) , {mj , %} 5 ma ,_ ma >m6 , _ mb . (12) 

Note that some coherence is returned to the ground state as a result of the "in terms". The 
derivation of the terms in the right-hand side of Eq. (10) is done by tracing over the field 
states with one-photon occupation number, a procedure which has been used in the case of 
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single multilevel atoms [see for example Ref. [16]]. The first two terms describe repopulation 
and recoherence of the ground manifold from the excited state manifold in a single atom. 
Cross coupling between atoms is reflected in the next two terms. Using again the solutions 
of Eqs. (7), the right-hand side of the in-term equations can be expressed in terms of ground 
state density matrix elements. A complete system of 16 linear equations is thus obtained by 
adding the in-term to the out-term contributions. 



III. SMALL SEPARATION (R 21 < A) 

In this limit, owing to angular momentum conservation rules, the propagators Gij vanish 
except for % = j. A few photon exchange processes between atoms are illustrated in Fig. 2(c), 
along with their accompanying propagators. Notice that due to momentum conservation the 
polarization of the emitted photon matches the polarization of the absorbed photon. Taking 
as an example the transfer of excitation from atom 1 in state |a) (with atom 2 in state ||)) 
to atom 2 in state |a) (with atom 1 in state ||)) depicted in Fig. 2, from Eq. (10), one finds 
that the propagator 

P — a r P 

< a ) | .... r^...* 

l — — T l —i-T 

p — -~r a r p— — a 

(b) | ....^OO.^ ; 

p — -r- a r p — -^r— a 

(c) ^n. i> 

i— — T l-^- T 



FIG. 2: Illustration of a few possible photon exchanges between atoms when the interatomic separation is 
small. Owing to angular momentum conservation, only diagonal elements of G are present. In (a) and (b) 
a transfer of excitation involving Am = and, in (c), Am = 1, transitions is shown. 
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A. Coupled basis 



As in [7], the two- atom system can be described by superpositions of states that are 
either symmetrical or antisymmetrical under particle exchange (Dicke states). The indis- 
tinguishability of the particles restricts the system to the symmetric subspace. The ground 
state manifold is symmetrized as follows: 



Ifi > = ITT), 

lflb> = ^(lU> + UT», 
|ff-i) = |U), 



(13) 



while ground-excited symmetric states are defined as: 

|ri> = -^(|oT> + |To», 

|r_i) = -^(|/?T> + lT0», 

|si> = ^(l«l) + ll«)), 

> = (|/U> + U £>) ■ 
Rewriting Eqs. (7) in terms of the new coefficients ri ; _i and si i , one finds 



(14) 



S-l 



~7Si -iAst ~\l r -i + ^X9-i, 
4 



(15) 



with quasistatic solutions 



5 -A 1 

-7r_! -«Ar_! - -7S1 , 



si 



ri 



3y/2 X (f 7 + *A) 
(|7 + ^A) 2 -I 7 : 

lA + f 7 



5-i , 



(16) 



r_i =0 
s_i = 0. 



XV 

a; 
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Notice that the symmetric states containing the excited (5 states have either identically zero 
or negligibly small amplitude, which substantially simplifies the calculation. 

In the following two subsections, these expressions for the state amplitudes are used to 
derive the time evolution of the collective coherence and population. In the new basis the 
matrix elements are denoted by , with i,j = —1, 0, 1. 



1. Coherence decay 

The expression for the collective coherence operator expectation value in the new basis 
is given by 

(P iT ) = V2(p 10 +po,-i )• (17) 
In the limit 7 < A, the following rate equations are obtained for the out-terms [from Eq. 
(9)] 

.out 4 x 2 

P10 = -g7^Pio , (18) 



■ out 14 x 2 
T 7 A" 2 



Po,-i = ~— 7t^Po,-i , 



while the in-terms, obtained from Eq. (10), evolve as 

4 X 2 



Pw = oTTaA) -1 > ( 19 ) 



3 A 2 

,2 



P0-1 = o^T^Po.-i 



X_ 

3 ' A 2 

Adding the out-term contribution to the in-terms, and with the notation 7 op = 7^ 2 /A 2 (op- 
tical pumping rate), the equations for the density matrix elements relevant for the coherence 
decay are 

4 

P10 = -g7o P (P10 - Po -1 ) , (20) 
P0-1 = -^lopPo-i ■ 
Substituting the solutions of the Eqs. (20) into Eq. (17), one obtains 

(At (*)> = l^ 2lopt (-1 + ^ hopt ) ■ (21) 

This is to be compared with the independent atom coherence decay derived in Appendix B 
[Eq. B4] 

(At (t)) md = e~ 7 < 
10 



The increase in the coherence decay rate for intermediate times, as shown in Fig. 3, can be 



coupled system 
■ independent atoms 



y t 



FIG. 3: Coupled system coherence decay for i?2l A vs. independent atoms decoherence. 



understood in terms of the exchange processes in the uncoupled basis illustrated in Fig. 2. 
For an independent atom, the state \a) is reached from the ||) state through the action of 
the classical field, and decays into ||) and |f) with rates 2(27/3) and 2(7/3), respectively. 
In the cooperative case, other channels responsible for coherence generation or decay appear 
owing to the presence of the second atom. 

An interesting behavior is observed at the initiation stage of the decoherence process, 
where the coupled system decoheres at a rate equal to that for independent atoms. However, 
this is not a general result, but rather a consequence of the initial state prepared with equal 
ground substate populations. When both atoms start in the same arbitrary state a \ l) + b |f) 
with a 7^ b and a 2 + b 2 = 1, the evolution of the coupled system coherence takes the following 
form 

(22) 



(P iT (t)) = 2abe- 2 ~ 1 °v t -a 2 + (1 + a 2 )^^ 1 
while the independent atom coherence evolves as: 

(P iT [t)) md = 2abe-^. 

Expanding the exponentials in Eqs. (22) and (23) for small times 7 op t C 1, we find 

(P iT (t)) md ^2ab(l- lop t), 



(23) 



(At (*)> - 2a6 



1 + 



b 2 -a 2 



lopt 



(24) 
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which shows that the decoherence rate of the coupled system is modified by the term (b 2 — 
a 2 )/3. This indicates that, given a population imbalance between the up and down states, 
at the initiation stage, the decoherence rate of the coupled system can be either larger or 
smaller than the one for the independent atoms, and vanishes for the balanced case only, 
when a = b = 1/ y/2. 




FIG. 4: An " imbalanced" superposition state with coefficients a,b = 1/10, \/99/10 decays exponentially 
at a rate 7 op , under the independent atoms assumption. However, when the decay is cooperative, an 
imbalanced state with a < b decoheres at a faster rate than one with a > b. 



2. Population Transfer 

In the symmetric basis, the expectation value of the collective population operator is 
expressed as 

(P n )= Pu + ipoo • (25) 

Three density matrix elements are coupled to each other: pu , poo and . The evolution 

resulting from the classical field can be obtained [from Eq. (9)] as 

. out 

Pn =0, (26) 

. out 8 
POO = --^1 op POO , 

. out 20 

P-1,-1 = g-7opP-l -1 , 
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while the in-terms are given by [see Eq. (10)] 



. in 4 

3 

. in 4 



Pll = g7o P P00 



Poo = gTopP-i-i + ^loppoo , 

■ in 16 

P-i,-i = y7opP-i -l • 

Combining the in-terms with out-terms, one finds rate equations 

4 

Pll = g7o P P00 , 

Poo = g7o P (P-i ,-i - Poo ) , 
4 

P-i-i = -gTopP-i ,-i > 



(27) 



(28) 



that are solved to give 



(At (0) = i - \z~ hopt (3 + 2 7op i) ■ 



This is to be compared with the independent atom population evolution 



(P n (t)) md = 1 - e" 5^ f . 



(29) 



(30) 



where only one mechanism for populating state |f) is present: excitation of state \a) by the 
classical field followed by decay at a rate 7/3 to state 



Q_ 

V 



coupled system 
independent atoms 



FIG. 5: Cooperative optical pumping (from level |1) to level |f)) for i?2l "C A vs. independent atoms 
optical pumping. 
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The increase in the population transfer rate, plotted in Fig. 5 vs. the independent atoms 
rate, is due to photons emitted by the second atom that drive the first atom into the \a) 
state. A faster excitation of state \ct) leads to a faster transfer to |f) for intermediate times. 

B. Mechanisms for coherence generation 

It is interesting to isolate the dynamics of a single atom with the purpose of identifying 
the mechanisms that lead to the modification of its radiative properties due to the presence 
of a second atom in its vicinity. We analyze the rate of change of the expectation value 
of the one atom coherence operator = | X) x ( T 1 1 - This is expressed in terms of density 
matrix elements as (^ a ^^) = Ptt;lt + PlUll > anc ^ ^ ^ s f° un d to satisfy the following equation 
of motion [using Eqs. (9) and (10)] 



The first term in the right-hand side in the above equation simply indicates the decay of 
the coherence of the independent atom at the expected rate 7 op . The second term contains 
the modification induced by the action of the neighboring atom. At the moment when the 
interaction between atoms is initiated, the density matrix of the coupled system can be 
factorized and this term can be written as: —^op^i 1 ^ ^ci 2 ^, where = — 
is the population difference operator for the first atom. The significance of this term 
is that x polarization (coherence) established in the second atom induces x polarization in 
the first atom through the vacuum, given a population difference. 

We proceed now to analyze the origin of this coupling term by considering two distinct 
situations in which atom 2 is prepared in a superposition a | j) + b |f) exhibiting coherence 
equal to ab, while atom 1 is prepared either in the |f) state or the ||) state. In both cases 
no initial x atomic polarization in the atom of interest is present. When starting with the 
population in the |f) state, using Eqs. (9) and (10), expressions for the density matrix 
elements present in the right-hand side of Eq. (31) can be derived, and an out-term is found 
to be responsible with the generation of coherence 




(31) 




(32) 




d 
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The process leading to this can be represented as follows 

X 7 X / Q Q\ 

PfTiU *■ PTT;T« > Ptt;«t > Ptt;lt > l dd J 

where in the first step, the field induces an excitation in the second atom, followed by a 
swap of excitation between atoms through the vacuum (collision-like effect) and a stimulated 
emission from the first atom. In the second case, where the population is initially stored in 
the ||) state, both an in-term and an out-term are present: 

Jt^^T = -\^ (34) 

di\ a - {t) ) = r° pab - 

The coherence generation through the out-term is similar to the process shown above [Eq. 
(33)], while the in-term takes the following path 

where consecutive excitations for both atoms are followed by spontaneous decay into a state 
with coherence in the first atom. 

IV. ARBITRARY SEPARATION 

Simple analytical results are not available in this regime. Coupling through propagators 
other than Gu takes place. The polarization of the emitted photon doesn't have to match 
the one of the absorbed photon, a situation which is illustrated in Fig. 6, where the coupling 
of the two- atom state \a f) to other states is shown and the corresponding elements of the 
G matrix (including non-diagonal ones) responsible for the coupling specified. 

A. Numerical Results 

The calculations are now performed in the uncoupled basis. Numerical solutions of 12 
coupled rate equations give the coherence decay, whereas 13 rate equations are solved to 
obtain the population transfer rate. The output of our numerical simulations is dependent 
both on time and on the spherical coordinates of the second atom Ri\,Q and <p. Since 
the system has azimuthal symmetry, the interesting cases are obtained by varying i? 2 i and 
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FIG. 6: The transfer of excitation stored in the coupled state \a |) can take place through a channel 
governed by a diagonal propagator Goo represented in ( a )j °r through non-diagonal propagators Goi ,Gio and 
Gi-i as in (b), (c) and (d). 

9. Two orientations of the system with respect to the field propagation are examined: 
R 21 || k L (0 = 0) and R 21 1 k L (6 = tt/2). 

1. Coherence Decay 



The 4 density matrix elements responsible for the coherence [Eq. (6)], 
are coupled to 8 more through nondiagonal elements of G. A sys- 

tem of 12 linear differential equations has to be solved containing 

{PlUll ;PU;U >ALT;U >All;T4. jAUUT > PlV,U 'PU;U ; PtUIT > PlT;lT >ALT;U >PTT;U >PTT;It}- Tne 
expressions for all 12 rate equations are not given here; instead a single one is written to 

illustrate the way the coupling among different states comes into play 



lop 

3 



~ 5 PU;U 



+e 



-ikx J R2i cos ( 



V2G 10 (R 2 i,9,ip)pu-il +2G 00 (R 2 i,0,(p)p l1 
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+2G n (R 21 ,6,<p)p mi 

,v, l/) ., j( ,,s« j 2Gl1 ( i? 2i,^,^)pn;U + V / 2G i ( J R 2 i,^,<^)pn; 

+V2G 01 (R 2 i,9,ip)p n . n 



(36)" 



The time evolution of the collective coherence is shown in Fig. 7 for kiR 2 i = 0.7. For 
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FIG. 7: The two-atom coherence for fe^i^l = 0.7 is plotted as a function of time for independent atoms, 
R-21 || ki and R.21 -L k^, respectively. 

both R21 || and R21 _L kx,, the coherence decays more rapidly than for the independent 
atom case, and the decay for R 2 i || k L is faster than that for R 21 _L k L . To obtain some 
idea of the dependence of the decay rate on interatomic separation, we plot in Fig. 8 the 
decoherence rate (( P^(t) ) / (P^(t))) as a function of distance, at a fixed time t = l/7 op , 



for both the R 2i || k L and R 2 i _L k L . Even though a numerical solution has been used 
to obtain the plot in Fig. 8, an approximate analytical treatment can provide insight into 
the qualitative nature of the results. In particular it can help explain why the parallel case 
decay rate is larger than that for closely separated atoms. The coupled equations of motion 
for the coherence operators associated with each atom [similar to Eq. (31) but with the 
difference that now the coupling coefficients are i? 2 i dependent] are solved approximately at 
a fixed time. 

It is found that the perpendicular case differs from the close atoms case only insofar as the 
coupling between atoms is modulated by the real part of G 00 (R 2 i , tt/2,0). An analytical 
expression for a fixed time t = l/jop gives a decoherence rate of the collective coherence 
that varies with the separation as 7 op [l + 0.21 Re(G o (-R21, tt/2, 0))]. 

There is, however, a fundamental difference between the perpendicular and parallel case 
that is reflected in the system's response. In the perpendicular case the spatial phase of 
the laser field does not enter since k^ • R 2i = 0. As a result the spin coherence associated 
with each atom evolves in an identical fashion. On the other hand, in the parallel case there 
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is a relative phase difference of k L • R 2 i associated with the laser field as it interacts with 
atoms at the two sites. Consequently, the response of the two atoms need no longer be 
identical, since the spatial symmetry has been broken by the field. Three coupling terms are 
present in the parallel case: Re(G o (R21, 0, 0) cos /CL-R21), Re(G o (R21, 0, 0) sin k^R^) and 
Re(Gn (R 2 i, 0, 0) sin/cii?2i); the importance of the laser induced spatial phase is evident in 
these expressions. Owing to the extra couplings, the two atoms accumulate different spatial 
phases, and the collective coherence (obtained as the sum of individual coherences) shows a 
spatial modulation that varies as Re(Gn (R21, 0, 0) sin k L R 2 i). A full analytical solution for 
the variation of the decay rate with the distance for any fixed time is not available; however, 
a perturbative treatment for small times {^ pt <C 1) indicates an increase in the decay rate 
2/9 Re(G 2 n (R 21 , 0, 0) sin 2 k L R 21 h 2 op t. 




FIG. 8: For a fixed time t = the decoherence rate of the two atoms (in units of 7 op ) is plotted for 

separations varying from to R 2 ± = 10/ k^. The dotted line respresents the independent atom decoherence 
rate (j op ). 

Just as in the case of close atoms, we extend our simulations to analyze the decay of an 
arbitrary initial state a \ l) + b The observed behavior is quite different here. Even for 
relatively large separation (i? 2 i = A/27r), a state with most population in the down state 
decays much faster than in the independent atom case, while an inhibition of decoherence is 
obtained when the initial state is prepared with more population in the up state [as shown 
in Fig. 9 for R 2 i || k L and R 2i _L k L ]. 
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FIG. 9: The decay of the collective coherence of an " imbalanced" superposition state with coefficients 
1/10 and \/99/10 is plotted above. Enhancement of the decay rate can be obtained when the initial state 
has more population in the up state, while the opposite case gives rise to an inhibition of decoherence. 

2. Optical Pumping 

One more density matrix element (ptT;TT ) ls coupled to the 12 listed before and a system 
of 13 rate equations is solved to obtain the collective population transfer as a function of 
time. 
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FIG. 10: The population transfer for k^^X = 0.7 is plotted for independent atoms, R21 || and 
R21 _L k^, respectively. 

An interesting behavior is obtained in both the parallel and perpendicular cases (see Fig. 
10): the optical pumping rate is initially larger and afterwards smaller than the one for 
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independent atoms. The saturation effect at large times is due to the addition of a decay 
channel (the \(3) state) which provides a way for the transfer of population from |f) back to 
the initial ||) state. This accounts for a slow down at large times where the repopulation 
from the second atom (resulting in population in state \/3)) is comparable with the population 
produced by the classical field. 

V. POLARIZATION SWAP 

As seen in Sec. Ill B, a partial transfer of coherence from an x polarized atom to an 
initially unpolarized atom (in a z state, either up or down) can be achieved. The two distinct 
situations discussed there can be extended for variable interatomic separations. Figures 11 
and 12 show the evolution of the coherence of the first atom as a function of time for 
R21 <C A and R 2 i = \/2n and \/ir, respectively, when R 2 i || k L . In both cases, owing to the 
oscillating nature of the coupling between atoms (as a function of separation) the sign of 
the effect produced in the initially z polarized atom depends dramatically on the distance. 
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FIG. 11: An x polarized atom 2 induces an x polarization in an initially z polarized atom 1 (spin down). 
The coherence in atom 1 is plotted as a function of time at R21 <C A, at R21 = and at R21 = 2/kz,, 

for the case when R21 || k^. 

The decay dynamics are still given by Eqs. (32) and (34), with the distinction that the 
decay parameter, which for close atoms is simply equal to 7 op /3, is now spatially modulated 
by Coo(-R2i, 0, 0) cos(/cii?2i)- In the case where the first atom is initially in the up state, 
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the coherence is driven only by the out-term, which leads to a simple time behavior, where 
(Pjf(t)) changes sign only due to the spatial oscillation of the coupling term. In contrast, in 
the down case, the competition between the in-term and out-term leads to a change in the 
sign of (Pjf(t)) for intermediate distances. 
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FIG. 12: Atom 1 starts here in state up. The polarization transfer is plotted as a function of time at 
R21 <C A, at R21 = 1/Jcl and R21 = 2/kL, for the case when R21 || III- 



VI. CONCLUSIONS 



We examined a system of two multilevel atoms coherently driven by a single-mode classi- 
cal laser field and coupled to the electromagnetic vacuum. The decoherence (of a quantum 
superposition stored in the ground sublevels) that necessarily accompanies the process of 
manipulation of the atomic states has been analyzed in the context of cooperative behavior. 
Two cases have been treated, where the atoms are either at the same position or separated 
by a distance comparable to the optical wavelength. It has been found (not surprisingly) 
that, for the case of close atoms, the "communication" between atoms leads to an increase 
in both the decoherence and population transfer rates. With increasing interatomic separa- 
tion, in the case of the field propagating perpendicularly to the line joining the atoms, the 
decoherence rate is less than that for close atoms. This is an expected result since, from 
the point of view of the classical field the atoms are located at equivalent positions, and a 
simple decrease of the interatomic coupling due to the increasing separation is expected. For 
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certain distances, however, owing to the oscillating behavior of this coupling, a small effect 
of decoherence inhibition is also observed. A more interesting situation arises for the case 
when the atoms are aligned parallel to the field propagation direction. The equivalency of 
positions does not hold here anymore, and a spatial phase difference between atoms resulting 
from the classical field is established. The coupling through the vacuum is modulated by this 
spatial phase difference, and a considerable enhancement in the decay rate is observed at 
separations of order \/2n. These results will be generalized in a future planned publication 
to a large ensemble of atoms in a pencil-shaped geometry. In such a medium, at Fresnel 
numbers close to unity, the atoms are practically aligned along the direction of the field; for 
large optical depths, the phase effect described above is expected to lead to a substantial 
increase in the decay of the collective atomic coherence. 
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VIII. APPENDIX A: EXPLICIT EXPRESSIONS FOR PROPAGATORS 

The expressions for the propagators involve spherical harmonics and Hankel functions of 
the first kind, h (k R) and h 2 {k R) [Ref. [14]]: 




Gn(R) = V^h (k R)Y ,o(R) - -J—h 2 (k R)Y 2fi (R), 

A V 

Goo(R) = V^ho(koR)Y 0fi (R) + ^h 2 (k R)Y 2fi (R), 



(Al) 



G 



-i-i 



(R) = G lil (R); G 0) -i(R) 



Gi,o(R); Go,i(R) 



G_ 1)0 (R). 
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IX. APPENDIX B: DERIVATION OF THE DECOHERENCE AND POPULA- 
TION TRANSFER RATES FOR A SINGLE 4-LEVEL ATOM 



We explicitly derive the equations of motion for the ground state coherence and popula- 
tions for a single 4-level atom interacting with a a + polarized field. In the perturbative limit 
where a maximum of one excitation is allowed, a basis set in the Hilbert space of atom and 
radiation field is comprised of states ||) |0), |T) |0) , |a) |0), \f3) |0), ||) |k, A) and |t) |k, A). 
The first ket denotes the state of the atom while, the second one describes a vacuum field 
or a one photon state with wave number k and polarization A. The state amplitudes obey 
the following equations of motion: 

K-fi = -lK-,0 -i&b a . +ixb i0 , (Bl) 

&/3;0 = -7&/3;0 ~ lA&fto , 

bv,o = 0, 

bl;k,X = igj^x ba;0 , 

Using a master equation approach one can write density matrix equations of motion for the 
ground state sublevels as: 




The first observation that we make here is that the terms in the right hand side of the above 
equations, that are due to the field (out-terms) and to the coupling to the vacuum (in-terms) 
can be derived separately. In other words, the presence of dissipation can be neglected when 
writing equations describing the driving effect of the field and later added to the equations 
phenomenologically. The derivation of rate equations can now be carried out by writing 
equations for the ground-excited coherences and excited state populations and coherences 
and adiabatically eliminating them. However, this is an unnecessary complication; the 
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second observation that we make provides an easier way of doing this. Only the first 4 
equations in Eqs. (Bl) have to be solved and the excited amplitudes can be written in 
terms of ground state amplitudes. Next, the density matrix elements can be written simply 
as if the state of the system were a pure state (as products of amplitudes) and replaced 
in Eqs. (2) to obtain rate equations. The separation of in-terms from out-terms in these 
equations insures the validity of this approach. Following this recipe, it is found that (in the 
limit A > 7) 

(B3) 
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Replacing these expressions into Eqs. (B2), rate equations are obtained in a final form: 

2 

Pll = gTapPu ' ( B4 ) 

Pll = -lopPll + l-^PU , 

4 2 
Pa = -fi°pPn - 2 ^o P pii = -gTapPu • 

The decoherence of an initial state 775 (||) + |T)) can now De calculated as described by the 
evolution of (neglecting the phase associated with the AC Stark shift) 

PIT (*) = \r"»\ (B5) 
while the population transfer from state ||) to state | j) is given by 

p n (t) = l-e-h^. (B6) 
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